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I. INTRODUCTION 

g-state Potts models have played an important role in 
condensed matter physics,!]. Here we study geometrical 
aspects of the critical Potts clusters, in two dimensions. 
The g-state Potts model is defined through the Hamil- 
tonian 

w = -A'E('5-'^-. -1)' (1) 

where denotes the summation over nearest neigh- 

bor sites the spin variable cr^ can take any of the 
values 1, 2, . . . , g and K is the thermal coupling with the 
factor l/k-gT absorbed in it. It is possible to extend q 
to real values 0, but here we concentrate on the Potts 
model with integer values of q. In particular, we study 
Potts models with q = 1,2,3, and 4, where the ther- 
mal phase transition at the critical inverse temperature 
Kc = In (1 -I- ^/g) P] is of second order. Below we present 
simulations at K = Kc- 

One defines the fractal clusters in the Potts model 
through the Fortuin-Kasteleyn (FK) Q cluster decom- 
position, which states that the model can be mapped 
onto a general percolation model. The partition function 
of the Potts model Z = Tr^-e^ can be expressed in terms 
of bond variables as Z = Trbondsp''(l ~ py'q^", where b 
is the number of bonds and n is the number of interac- 
tions that did not form a bond in a configuration with 
Nc clusters Here, p = 1 — e^^, and Trbonds means 
a summation over bonds. Thus, the problem of a ther- 
mal lattice model can be mapped to a graph problem. 
The FK decomposition has been the starting point for 
efficient cluster algorithms ^ Q for simulation of spin 
models. 

The Potts model has been shown to exhibit a rich crit- 
ical behavior, and it has been related to a number of 
problems in lattice statistics 0]. Although of great the- 
oretical interest in itself, it also has many experimental 
realizations. The 1-state Potts model is equivalent to 
a bond percolation problem and the 2-state Potts 
model is the same as the Ising model 0. The q = 3 



Potts model has been shown to describe absorbed mono- 
layers on two-dimensional (2D) lattices 7, 8J. Domany 
et al. suggested that N2 absorbed on krypton-plated 
graphite should exhibit the same critical behavior as the 
q — 4: Potts model. More references of the experimental 
realizations can be found in the review article by Wu . 

Here we study geometrical aspects of the critical Potts 
clusters in two dimensions. This is in direct analogy 
with the geometry of percolation clusters, which has been 
widely studied [Tfl lO, U2, ■ Specifically, we measure 
the fractal dimensions Dm, Dh, Dep,Dsc, and Dq 
describing the scaling of the cluster's mass, hull, external 
accessible perimeter jl^ , singly connected bonds 0| and 
the gates to narrow-gate fjords |0, respectively, with its 
radius of gyration R. As emphasized by Coniglio [l5l| . 
many of these fractal dimensions have been derived ana- 
lyti cally pdl | . Some others have been found more recently 
[12. Ha l . Although there exists much numerical work on 
the percolation clusters (i.e. g = 1), we are not aware 
of any detailed numerical study of most of the above 
mentioned quantities for g > 1, and especially when q 
approaches the critical value qc — 4. 

Section^]describes the numerical methods used in the 
simulation of the Potts models. Our numerical simula- 
tions show that the asymptotic power law dependence of 
the various masses on R is approached relatively slowly, 
and therefore the analysis of the data must include cor- 
rection terms, particularly as q approaches qc- The the- 
ory developed to obtain these correction terms [l^l is 
briefly summarized in Section IIIII We compare our nu- 
merical data with the exact predictions in Section IIVI 
Finally, we present the summary and conclusions in Sec- 
tion |V| 



II. SIMULATION 

Numerical simulations of spin models have developed 
from the local spin flip type algorithm to the more 
advanced cluster algorithms 0, Q ■ Our simulations were 
done on a 2D square lattice with both open and periodic 
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boundary conditions. Clusters of the g-state Potts model 
were generated using the Swendsen-Wang algorithm 0, 
which is based on the cluster decomposition by Fortuin 
and Kasteleyn [sj. The size of the system in our simu- 
lations was 4096^ spins for all q. Figure shows sample 
clusters for different values of q. 

All simulations were started with a homogeneous ini- 
tial condition, with all spins initially parallel to each 
other. First we thermalized the system to allow the 
model to equilibrate. Thermalization was checked by 
measuring both the energy per spin e, directly from the 
Potts Hamiltonian, and the magnetization per spin m, 
using the representation of Potts spins in a, q ~ 1 dimen- 
sional space p. 

Thermalization of large spin systems takes a very long 
time. The quantities of interest in this work show a rel- 
atively slow approach to the asymptotic values. Thus, 
extremely large lattices are required for the analysis 
of the scaling behavior of the cluster subset masses. 
When performing simulations on lattices of linear size 
L = = 4096, about 20000 Monte Carlo steps (lattice 
sweeps) are needed to equilibrate the system. 

We devised a simple method to overcome the problem 
with long thermalization times. We started the thermal- 
ization with a small lattice of size Li and thermalized 
it. We then periodically copied the spin configuration of 
the small lattice to a lattice with twice as large a size, 
L2 — 2 X Li, and thermalized it. We continued this 
process until the desired system size was reached. In 
practice, it is recommended to compare the values of e 
and m obtained this way with the values obtained from 
conventional thermalization to be sure that the system is 
really thermalized. Alternatively, one can continue run- 
ning the simulation and collect the values of thcrmody- 
namical variables as a function of time and check that 
there is no increasing or decreasing trend in them. The 
thermalization method described above allows a speed- 
up by an order of magnitude in thermalization for a Potts 
spin system of 1024^ spins. 

After the spin system was thermalized, we took sam- 
ples of the cluster configurations after every 20 spin up- 
dates (corresponding roughly to the correlation time for 
the present system sizes) to get uncorrelated samples. 
Each cluster of the present configuration was taken sep- 
arately under investigation. However, as a precaution to 
avoid some of the finite size effects, we collected the data 
only from those clusters which did not involve spins on 
the boundary. For the remaining clusters we determined 
the masses of the cluster subsets. The total mass of the 
cluster is defined as the number of occupied bonds in 
the cluster. The number of bonds belonging to the hull, 
external perimeter and the number of singly connected 
bonds were counted using directed walkers that walk on 
the appropriate cluster perimeter 0, ^ . 

To measure the number of gates to fjords of different 
gate sizes (Sg) a walk was initiated on the EP. The walker 
starts from the left vacant neighbor of the lower-left site 
belonging to the cluster and goes around the cluster on 



the EP always trying to turn to the right. At each site, 
we look at the neighbors on the left hand side and check 
whether the sites within a predefined distance Sq belong 
to the EP or not. If the sites (and bonds) up to the 
distance Sq are not EP sites and the site at the distance 
Sg belongs to the EP, the walker is about to enter a 
fjord with a gate size Sq- All fjords with different gate 
sizes were counted during a single walk around the cluster 
with an array of boolean variables that indicate whether 
the walker is within a fjord of given size. If the walker 
was in a fjord of a gate size Sg, fjords with gate size 
S'g > Sg were not allowed in the statistics. In practice, 
all measured values of Sg gave similar results, and we 
report only the results for S'g = 1. 

As noted by Aizenman et al. 12^ , the scaling concerns 
fjords whose size Lp is comparable to that of the clus- 
ter, R, and whose gate's width is much smaller than Lp- 
The reason for this is easy to understand, since small 
kinks and pits are a natural part of the fractal cluster's 
perimeter. Thus, only large enough fjords were included 
in the statistics. This was taken into account by choos- 
ing a suitable parameter s, and counting only fjords with 
Lp > sR. 

Since any thermalized spin configuration contains clus- 
ters of many sizes, we collected the data in multiplica- 
tively increasing bins of the size of the cluster. Each bin 
contained the clusters of sizes within [Ri,Ri+i = \f2 

III. CORRECTIONS TO SCALING 

In our recent publication 'iT"] , we derived theoretically 
the corrections-to-scaling terms for the various cluster 
subset masses. The corrections arise from three different 
sources: (a) using the renormalization group approach, 
(b) mapping to the coulomb gas, and (c) considering the 
uncertainty of the correct measure for the linear cluster 
size, which implies corrections of the order of 1/R. 

The first correction relates to the dilution field ip, 
which is generated under renormalization even when one 
starts with the non-diluted case 0. Solving the renor- 
malization group recursion relations (RGRR's) for ip{£), 
where £ is the scale variable, substituting the solution to 
the RGRR's for the field hs conjungate to the density 
ps — Ms/R'^ and finally calculating the scaling of the 
density ps{£) with the cluster's linear size R = yields 
the following predictions for the approach of each mass 
Ms to the asymptotics. Here, S = M, H, EP, SC or 
G. 

In the q — A case, the renormalization group calcu- 
lation is exact, yielding logarithmic corrections to the 
scaling of Ms{R): 

Ms oc R°''(\ogR + BlogilogR) + E)-''''/'' 

x(l + 0(loglogi?/logi?)), (2) 

with Ds = ys{q — 4) and cs/a as given in Table |I| (see 
below) . Note that B is universal, and the non-universal 
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constant E is the same for all S. Equation (O generalizes 
the logarithmic corrections of Cardy et al. jl9i| . 

We were recently directed to a previous study by Van- 
derzande and Marko "2^, in which they considered cor- 
rections to scaling for pcrcolativc properties of the <? = 4 
Potts model. They also found logarithmic corrections to 
scaling in aggreement with those presented here. 

For g < 4, to leading order in e' = ^4 — q, the same 
procedure yields 

Ms oc - Si?-'')-"-'/" « + /si?"^), (3) 

where Ds ^ ys ^ cse' and 6 « 2ae'. Note that to the 
lowest order in e', the ratios fs/fs' are universal, being 
equal to cs/cs'- This is similar to analogous ratios for 
thermodynamic properties in the usual e-expansion j2ll |. 
This universality should hold to all orders in e' . Expand- 
ing the exact Ds (Tabled in e' yields cg. Using also 
a = I/tt [t^IT^ yields our predictions for cs/a (given in 
Table P), to be used in the fitting procedure. 

The second source of corrections involves new contri- 
butions to the relevant pair correlation functions in the 
Coulomb gas representations, den Nijs "22] derived such 
corrections to the order parameter correlation function. 
Since correlation exponents x are related to the fractal 
dimension via D = d — x, the correction exponents can 
be related to the corresponding correction terms for the 
scaling of the mass Mm, yielding the leading correction 

0' = 4/5, (4) 

where g is the (g-dependent) Coulomb Gas coupling con- 
stant (see Table nj. 

Using a similar approach we found in the case of the 
hull and the singly connected bonds that the leading cor- 
rection exponent is given by 

0" = 2/g. (5) 

We argued that this correction would also hold for the 
external perimeter [T^ . 

The last source of corrections involves 'analytic' terms, 
coming e. g. from linear cuts with dimensions (Dg — 1), 
I^S*! or from replacing i? by (R + A), since there are many 
possible candidates for the correct linear measure of the 
cluster. These would imply corrections of relative size 
1/R. 

IV. RESULTS 

We now present the numerical data from large scale 
Monte Carlo simulations of the Potts models. Our aim 
is to confirm the exact predictions of the fractal dimen- 
sions Ds in the cases where they are available and to give 
numerical estimates for the exponents that have not yet 
been calculated exactly. In addition, we want to numeri- 
cally confirm the corrections-to-scaling theory presented 
in the previous section. 



We obtain good agreement with the theoretically pre- 
dicted values for most of the fractal dimensions Ds- The 
worst agreement is found for the exponent of the external 
perimeter Dep for q > 2 Potts models. The reasons for 
this will be discussed below. However, fixing the correc- 
tion terms and performing fits only to the amplitudes and 
to the fractal dimensions Ds in the logarithmic deriva- 
tives of Eqs. (O and yielded estimates for the subset 
fractal dimensions that agree to the precision of 0.05 or 
better with the theoretical predictions of Table 

We are not able to give a quantitative numerical proof 
of the values predicted for the correction terms. Quali- 
tatively, our numerically evaluated functions Ms{R) dis- 
play a complex corrections-to-scaling behavior that re- 
quires more than one correction term. The theoretical 
predictions can be fitted reasonably well to the data keep- 
ing the important quantities such as 6, 9' (or 9") and cs /a 
fixed while leaving the amplitudes of individual correc- 
tion terms free. 

We start this section by studying the fractal geometry 
of clusters in the site percolation model, which is com- 
putationally easier to simulate. The reasons behind the 
difficulties in the comparison of the numerical data with 
the analytical predictions are discussed. We then proceed 
to present our numerical data on the Potts clusters. In 
all the figures of this Section, whenever the error bars are 
not shown they are smaller than the size of the symbols. 



A. Site Percolation Clusters 

For g = 1 we simulated site percolation on a square 
lattice of size 24576^, using the Newman-Ziff cluster la- 
beling which is an improved version of the Hoshen- 
Kopelman algorithm 12^. Thus the linear lattice size 
was 6 times larger than in the simulation of Potts clus- 
ters with q> \. 

To get a feeling of what kind of problems are present 
in the fitting procedure when the theoretically predicted 
correction terms of Eqs. lj2Jl and Q are fitted to the 
numerical data on the cluster subset masses Ms{R), let 
us consider as an example the scaling of the number of 
singly connected bonds. Figure[21illustrates the scaling of 
MsdR) with the cluster size i? on a double logarithmic 
scale. The solid line in the main figure indicates the 
predicted slope. The fit to the data yields an estimate 
for the asymptotic fractal dimension Ds which is less 
than 0.01 off the exactly known value Dsc — 3/4. 

Although the asymptotic scaling regime Msc{R) oc 
R^^^ can be seen here, there are difficulties in the ex- 
traction of the correction terms of Eq. The small- 
est value of R included in the linear fit to the data on 
log — log scale in Fig. |21 corresponds to the regime where 
the influence of the correction terms are about to vanish, 
thus justifying fitting without any correction terms. The 
saturation to the asymptotics can be seen more clearly 
in the inset of Fig. [3 where data are scaled with the pre- 
dicted asymptotic behavior Msc{R) I R^^^ ■ The inset 
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shows that at about R ~ 300 the correction terms can be 
neglected in this case. However, at the same point the 
statistics becomes so noisy, making a precise estimation 
of the correction terms difficult. An additional difficulty 
arises from the fact that the finite size of the lattice is not 
taken into account in any way in the finite size scaling 
form of Eq. ©. Due to the finite system sizes, statis- 
tics of the large cluster is biased in such a way that only 
the compact clusters fit in the lattice without touching 
the boundaries. The extended clusters having for exam- 
ple more EP sites than compact clusters with the same 
radius of gyration R, are absent. This bias cannot be 
taken into account by any known correction terms. We 
tried to extrapolate the data from different system sizes 
to obtain an asymptotic curve for an infinitely large sys- 
tem, but the statistics is far from sufficient for such a 
procedure. 

B. Potts Clusters 

In the case of Potts clusters, system sizes that can be 
used in the simulations are much smaller than those in 
the site percolation case, since in addition to the spin 
variables, also bonds must be stored in the computer 
memory. This causes the finite size effects to be even 
more pronounced than those present in the site percola- 
tion model simulations. In addition, the correction ex- 
ponents in the 1 < g < 4 Potts models are smaller than 
those in the q = 1 case. Also, the logarithmic corrections 
present in the 4-state Potts model are weaker than any 
of the power law corrections in g < 4 models. Thus, the 
influence of the corrections-to-scaling terms extends to 
much larger values of R. 

The data analysis was done by fitting the theoretical 
predictions of Sect. Illll to the data. The nonlinear fit ting 
was done using the Levenberg-Marquardt method [26j. 
The measure for the quality of fits is x^- Values of 
close to one indicate a good fit. We determined the error 
bars of the fractal dimensions by fixing Ds to a range of 
values around the theoretically predicted one, and per- 
forming a fit to the amplitudes for each such value. The 
error bar on Ds was fixed as the range within which 
remained smaller than 2. One can perform the fit directly 
to Ms{R), or to the effective fractal dimensions D|^(i?) 
(which is the logarithmic derivative of Ms{R) [lol|V 

We found that for g < 4, fitting directly to the mass 
worked out better. If one wants to fit to the measured 
data directly, it is recommended to divide the measured 
data on the cluster subset masses by the exactly known 
asymptotic behavior to avoid problems with numerical 
accuracy. Thus, for g < 4, the fits are done to the form 

Ms/R''' = Es{l + fsR-' + fsR'"' + gs/R), (6) 

where Es, fsifs^9s,() and 9' are the fitting parameters, 
and the fractal dimensions Ds are kept fixed. 

In the q = A case, on the other hand, we found the 
fitting to the effective exponent Z??? better, and the fits 



were thus done to 

Df = dlogMs/dlogR 

^ ^ [csla)[B + C + \ogR) 

[C + log R){E + log R + B\og{C + log R)) 
+ZllogR. (7) 

Note that we have replaced the log log R term in Eq. Q 
by the more general log(C -t- logi?). Also, the logarith- 
mic derivative of the higher order term on the RHS of 
Eq. (0) was approximated by a simpler form 1/logi?. 
In the fitting procedure, the possibility of having many 
candidates for the correct linear measure of the cluster 
size was taken into account by allowing R to adjust to 
R + A (see Sect. Iin|) . 

Below, the results of our numerics are summarized for 
the various subsets. We do not go into the details of nu- 
merical estimates for the various amplitudes in the non- 
linear fits, since due to the relatively small range of the 
data and many fitting parameters the estimated error 
bars are large, and allow no comparison with, for exam- 
ple, the predicted amplitude ratios. Also, precise esti- 
mation of the correction exponents 0, 9' (or 9") as well 
as the parameters cs/a is impossible with the presently 
available range of data. Instead, we keep the correction 
exponents (or parameters) fixed and try to extrapolate 
the fractal dimensions Ds, and to demonstrate that the 
predicted forms of scaling in Eqs. {Tj) and (|SJ| are con- 
sistent with our numerical data. Specifically, in all the 
fits for 1 < (7 < 4, 6* and 9' (or 9") of Eq. ® were kept 
fixed (at the predicted values), while the amplitudes of 
each correction term were allowed to adjust. In the g = 4 
case, only cg/a was fixed in Eq. 10. 

The logarithmic corrections are most important for the 
singly connected bonds at g = 4, where theory predicts 
that Dsc = (see Tabled). Indeed, the solid line in 
Fig. O shows that a fit to Eq. {T)) is consistent with this 
theoretical prediction. In contrast, a fit with a single 
power law correction term 9 = 1/2 (dashed line in Fig. 
01 extrapolates to a wrong value near Dsc = 0.21! 



1. Mass 

Fig. ^ shows an example of the fit to the curve 
MM{q = 3)/R^'^' . In the fitting procedure, Dm, 9 and 
9' were kept fixed while Em, ImtI'm and qm were al- 
lowed to fit. The value x^ = 1.17 indicates that Eq. ® 
gives a good representation of the data. Our numerical 
estimates for the fractal dimensions Dm(q), determined 
as the range of values for which one has < 2, are 
1.90 ± 0.01, 1.87 ± 0.01, 1.85 ± 0.02 and 2.05 ± 0.15 for 
q = 1,2,3,4, respectively. These are in good agreement 
with the theoretical predictions. 
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2. Hull 

Figure|5lshows a fit to the number of the bonds belong- 
ing to the hull in the g = 4 Potts model. In this particular 
fit, Dh and Cn/a were kept fixed while B, C, Z and A 
were free to adjust. The value of = 1.11 indicates that 
Eq. {Tj) fits the data well. Our numerical estimates for 
the fractal dimensions Dniq) are 1.75±0.01, 1.66±0.01, 
1.59 ±0.03 and 1.50 ±0.01, for q 1,2,3,4, respectively. 
Agreement with the theoretical predictions is excellent as 
can be seen by comparison with the values in Table |l| 

3. External Perimeter 

Figure shows an example of the fit to the external 
perimeter data in the q — 2 Potts model. The fractal di- 
mension Dep and the correction exponents 9 and 0" were 
kept fixed at the predicted values while Eep, fsp, I'ep 
and gEP were free to adjust yielding = 1-77, implying 
a reasonably good agreement with Eq. Again, in 

the fits for q < 4, 9 and 9" of Eq. © were kept fixed 
and in the q = A case, cep/o, was fixed. The numer- 
ical estimates 1.33 ± 0.05, 1.36 ± 0.02, 1.40 ± 0.02 and 
1.48 ± 0.02 for q = 1,2, 3, 4, respectively, agree with the 
exact predictions. 

4- Singly connected bonds 

In Fig. Ul we show the number of singly connected 
bonds Msc{R)/ R^^^ against the cluster size R in the 
q = 2 Potts model. The fractal dimension Dsc and 
the correction exponents 9 and 9" were kept fixed to the 
predicted values while Esc, fscf'sc ^'^'^ Ssc were al- 
lowed to fit. The value = 1-21 implies good agreement 
with Eq. © . The numerical estimates for the fractal di- 
mensions are 0.75 ± 0.02, 0.55 ± 0.03, 0.35 ± 0.07, for 
q = 1,2, 3, 4, respectively. All the estimates for the frac- 
tal dimensions Dsc are in good agreement with the the- 
oretical predictions. However, the large value of « 3 
in the q — A case indicates some discrepancy between Eq. 
Q) and the data. 



scale was considered, since the scaling regime for the 
presently available cluster sizes is rather narrow. The 
parameter s governing the minimal ratio of the fjord 
size to the cluster size that we used was in the range 
0.1 < s < 0.2. The actual choice for the value of s does 
not affect the scaling law, but it merely determines the 
range where the power law behavior Mq ~ R^'^ starts 
(decreasing s shifts the maximum of the curves Mg{R) 
to the left). Our numerical estimates for |iI'G((7)| de- 
crease with increasing q. Our estimates for Dciq) are 
-0.90±0.05, -0.71±0.05, -0.63±0.05 and -0.59±0.05 
for q = 1,2,3,4, respectively. Our numerical estimates 
together with the theoretical predictions for all fractal 
dimensions Ds are summarized in Tabic Ull 



V. CONCLUSIONS 

The present paper examined the fractal geometry of 
the Potts clusters at the critical temperature. The aim 
was to find numerical evidence on the exactly derived 
subset fractal dimensions Ds [lllIT^ ITsllTill and to give 
estimates on the dimensions for which there is no exact 
prediction. We gave the first numerical estimate of the 
negative fractal dimensions Dq, describing the scaling of 
the gates to fjords [l^. 

Analysis of our numerical data revealed a slow and 
complex approach to the asymptotic behavior. If this 
is neglected in data analysis, wrong numerical estimates 
for the dimensions Ds follow. Using the corrections-to- 
scaling terms derived in our earlier publication in the 
fitting procedure, excellent agreement with most of the 
exact dimensions and data was found. The present qual- 
ity and range of data does not allow a unique quantitative 
confirmation of the exact correction parameters. 

To summarize, in the comparison between theory and 
numerics, extreme caution is needed in the extraction 
of the fractal dimensions Ds from the numerical data. 
The corrections-to-scaling theory presented already im- 
plies that the finite size effects arising from the finite 
cluster size are strong. In addition, effects coming from 
the finite lattice size lead to an uncontrollable bias that 
is very difficult to handle. 



5. Gates to Fjords 

Figure |H| shows our numerical data for the number of 
gates to narrow-gate fjords. The figure shows fits to the 
data along with the estimates for the fractal dimensions 
Dc- Our estimate Dc{q = 1) = —0.9 ± .05 agrees with 
the exact prediction Dc = -11/12 « -0.92 12]. Here, 
only a linear fit to the data on the double logarithmic 



VI. ACKNOWLEDGMENTS 

We would like to thank Carlo Vanderzande for drawing 
our attention to their earlier work. We acknowledge sup- 
port from the German-Israeli Foundation. This work has 
also been supported in part by the Academy of Finland 
through its Center of Excellence program. J. A. wishes 
to thank the Vaisala Foundation for financial support. 



[1] F. Y. Wu, Rev. Mod. Pliys. 54, 235 (1982). 



[2] R. B. Potts, Proc. Camb. Phil. Soc. 48, 106 (1952). 



6 



[3] C. M. Fortuin and P. W. Kasteleyn, Physica (Utrecht) 

57, 536 (1972). 
[4] R. Swendsen and J. Wang, Phys. Rev. Lett. 58, 86 

(1987). 

[5] U. Wolf, Phys. Rev. Lett. 62, 361 (1989). 
[6] E. Ising, Z. Phys. 21, 613 (1925). 
[7] S. Alexander, Phys. Lett. A 54, 353 (1975). 
[8] M. Bretz, Phys. Rev. Lett. 38, 501 (1977). 
[9] E. Domany, M. Schick, and J. S. Walker, Phys. Rev. Lett. 
38, 1148 (1977). 
[10] D. StaufFer and A. Aharony, Introduction to Percola- 
tion Theory, revised 2nd edition ed. (Taylor and Francis, 
Burgess Science Press, Basingstoke, London, 1994). 
[11] H. Saleur and B. Duplantier, Phys. Rev. Lett. 58, 2325 
(1987)). 

[12] M. Aizenman, B. Duplantier, and A. Aharony, Phys. Rev. 

Lett. 83, 1359 (1999). 
[13] T. Grossman and A. Aharony, J. Phys. A 19, L745 

(1986). 

[14] A. Coniglio, J. Phys. A 15, 3829 (1982). 

[15] A. Coniglio, Phys. Rev. Lett. 62, 3054 (1989). 

[16] B. Duplantier, Phys. Rev. Lett. 84, 1363 (2000). 

[17] A. Aharony and J. Asikainen, in Scaling and Disordered 
Systems, edited by F. Family, M. Daoud, H. J. Her- 
rmann and H. E. Stanley (World Scientific, Singapore, 
2002); will a lso app ear in Fractals (accepted for publish- 
ing) (2002); |cond-m at/0206367 

[18] N. Metropolis et al, J. Chem. Phys 21, 1087 (1953). 

[19] J. L. Cardy, M. Nauenberg, and D. J. Scalapino, Phys. 
Rev. B 22, 2560 (1980). 

[20] C. Vanderzande and J. F. Marko, J. Phys. A 26, 7391 
(1993). 

[21] A. Aharony and G. Ahlers, Phys. Rev. Lett. 44, 782 



(1980). 

[22] M. P. M. den Nijs, Phys. Rev. B 27, 1674 (1983). 

[23] Y. Gefen, B. B. Mandelbrot, A. Aharony and A. Kapit- 
ulnik, J. Stat. Phys. 36, 827 (1984); B. B. Mandelbrot, 
Y. Gefen, A. Aharony and J. Peyriere, J. Phys. A, 18, 
335 (1985) . 

[24] M. E. J. Newman and R. M. Ziff, Phys. Rev. Lett. 85, 
4104 (2000). 

[25] J. Hoshen and R. Kopelman, Phys. Rev. B 14, 3438 
(1976). 

[26] W. H. Press, S. A. Teukolsky, and W. T. V. B. P. Flan- 
nery. Numerical Recipies in C: The Art of Scientific 
Computing (Cambridge University Press, Great Britain, 
Cambridge, 1992). 



Ds 




g = 2 


q = 3 


g = 4 


cs/a 


a 


8 
3 


3 


10 
3 


4 




M (3 + 2)(5 + 6)/(8ff) 


91 
48 


15 

8 


28 
15 


15 

8 


1 

16 


H l + 2/g 
EP 1+5/8 


7 

4 
4 
3 


5 
3 
11 

8 


8 
5 

17 
12 


3 
2 
3 
2 


1 

4 

1 

4 


SC (35 + 4)(4-ff)/(8ff) 


3 

4 


13 

24 


7 

20 





-1 


e M,4-g)/9 


2 


4 
3 


4 
5 


(log) 




9' 4/g 
9" 2/g 


3 
2 
3 
4 


4 
3 
2 
3 


6 
5 
3 
5 


1 

1 

2 





TABLE I: Exact theoretical predictions. 
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S Ds{q = l) Ds{q = 2) Ds{q = S) Ds(g = 4) 



n e nenene 



M 


1.90(1) 


91 
48 


1.87(1) 


15 
8 


1.85(2) 


28 
15 


2.05(15) 


15 
8 


H 


1.75(1) 


7 

4 


1.66(1) 


5 
3 


1.59(3) 


8 
5 


1.50(1) 


3 
2 


EP 


1.33(5) 


4 
3 


1.36(2) 


11 

8 


1.40(2) 


17 

12 


1.48(2) 


3 
2 


SC 


0.75(2) 


3 
4 


0.55(3) 


13 

24 


0.35(7) 


7 
20 


0.03(8) 





G 


-0.90(5) 


11 
12 


-0.71(5) 




-0.63(5) 




-0.59(5) 





TABLE II: Comparison of the numerical estimates (n) for the 
subset fractal dimensions Ds with the exact predictions (e) 
where available. Uncertainties of the last decimal(s) for each 
Ds are given in parenthesis. 
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FIG. 1: Computer generated Potts clusters for (a) q = 1, (b) 
q — 2, (c) q = 3 and (d) q = 4 state Potts models. Colors in- 
dicate different subsets: SC bonds are shown in red, H bonds 
are shown in yellow and the rest of the bonds contributing to 
M are shown in blue. The EP bonds are colored green and the 
gates to fjords are marked by black circles, while the fjord is 
shown with a black line. For all the clusters, the total masses 
Mm are in the range 14400 — 17600. Note the decrease of Dh 
with q. 
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FIG. 2: Site percolation model. The number of the singly 

connected bonds Msc{R) vs. the cluster linear size R. The 
predicted slope Dsc ~ 3/4 is indicated by the solid line. 
The inset shows the scaled mass Msc{R)/R°'^'^- Note the 
saturation to the asymptotic scaling at i? « 300. 
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FIG. 3: Example of diflerent fits in the q — 4 Potts model 
to Msc{R)- The dashed line shows a fit with a single power 
law correction (6 = 1/2) and the solid line shows the fit to 
the logarithmic form of Eq. Q. Note the difference in the 
extrapolation to the R oo where the fits give Dsc = 0.21 ± 
0.01 with a single power law correction term whereas Dsc = 
0.03 ± 0.08 with the logarithmic form. 
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FIG. 4: Total mass of the cluster mass Mm in the 3-state 
Potts model. Solid line is the nonlinear fit to the data. For 
this particular fit ~ 1-17. 
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FIG. 5: Effective exponent -D^(-R) against 1/-R in the q = A 
Potts model. For the fit indicated by the sohd hne we find 
X' = 1-11- 




FIG. 6: Number of the external perimeter bonds, Mep versus 
1/R in the q = 2 Potts model. Solid line indicates the fit to 

the data (x^ = 1.77). 
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FIG. 7: Number of the singly connected bonds Msc against 
1/-R in the q = 2 Potts model. Solid line is the nonlinear fit 
for which = 1.21. 
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FIG. 8: Data for Mg{R), the number of gates to fjords on log- 
log scale. Different values of q arc represented by the symbols 
shown in the legends. Straight lines indicate the fits to the 
data; slopes give the exponents Da- 



